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1. Introduction 

The coherent manipulation of systems at their quantum level is an essential capability 
in many fields such as quantum information processing (QIP), quantum simulation and 
metrology. One aspect that separates a quantum system from its classical counterpart, is 
an exponential increase of state space as a function of the physical size of the system. The 
realization of arbitrary quantum operations on large systems is therefore a complicated 
task. Nevertheless, as in classical computation, any complex quantum operation can be 
approximately performed, with arbitrary accuracy, by concatenating operations taken 
from a finite set of simple gates. The finite set of gates from which any unitary operator 
can be constructed is called a universal gate set [1-3]. In the case of N two-level systems 
(a register of N-qubits), a universal gates set, for instance, may be comprised of a discrete 
and finite set of single qubit rotations and a two-qubit controlled-not (CNOT) gate [4], 
While it should be possible to apply each of these gates to any qubit, or qubit-pair, in 
the register, the number of gate types does not depend on the number of qubits involved. 

Trapped-ion qubits are one of the most promising platforms for QIP. Ions can be 
trapped for long times and laser-cooled to their ground motional state. Their internal 
electronic state as well as their motion can be controlled with high precision using 
lasers, microwave and radio-frequency fields [5,6]. Owing to the isolation of trapped 
ions from their surrounding, quantum states can have an extremely long coherence time 
[7]. A convenient choice of internal states for information encoding are the hyperfine 
and Zeeman states in the electronic ground manifold. These states are resilient against 
spontaneous decay and benefit from, commercially available, high-quality frequency 
sources. However, addressing individual ions and producing spin dependent forces 
presents great challenges due to the long wavelength involved. The common solution 
to this problem involves the use of Raman transitions, which implies relatively high 
intensity laser light. Another, recently proposed and demonstrated solution, is the use 
of near field microwave gradients [8]. A different choice of qubits are optical qubits, 
in which the two qubit levels are separated by an optical transition. Several of the 
heavier ions, that are used for QIP purposes, have low lying D levels (lower in energy 
than the P manifold) with a life-time of the order of one second due to the forbidden 
dipole transition to the S ground state. Here, the qubit levels are encoded in the S and 
D manifolds which can be coupled by a narrow optical electric-quadrupole transition. 
The transfer of the transition photon momentum to optical qubits is straightforward 
and so is individual ion-qubit addressing with spatial resolution that is limited by the 
optical radiation wavelength. These advantages are compensated by the fact that the 
coherence time of such qubits is limited by the linewidth of the laser that serves as a 
local oscillator. 

In this paper we present a complete set of operations on a two-qubit register realized 
of two trapped-ion optical qubits. These includes single qubit rotations and a two-qubit 
entangling gate. Our main focus here is on the quadrupole transition at 674 nm. A 
description of the ion trap and laser systems is detailed in [9]. The relevant internal states 
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Figure 1 . a) Energy level diagram and relevant transitions in 88 Sr + ion. Straight 
arrows indicate transitions that are driven by lasers, wiggled arrows indicate transitions 
due to spontaneous emission, where the thick wiggled arrows represents fluorescence 
which is used for detection, b) An example of a measured two-qubit fluorescence 
histogram. 


and transitions in a 88 Sr + ion are depicted in Fig.l. We encoded quantum information 
in two optically separated manifolds in the | S) = |<S'i/ 2 , m =+i/ 2 ) and | D) = |i? 5 / 2 ,m=+ 3 / 2 ) 
states. The two manifolds are coupled by a quadrupole transition where the excited D 5 / 2 
level has a lifetime of 390 ms. Driving coherent operations between the two manifolds 
was accomplished using a narrow linewidth 674 nm laser in a master-slave configuration. 
Due to the finite Lamb-Dick parameter 0.04 the motional degrees of freedoms 

of the ion can be accessed as well. This is a necessary capability since getting the two 
qubits to interact can be efficiently done only through their Coulomb interaction. The 
scheme that was implemented here for entangling two ions, is the Molmer-Sorcnscn (MS) 
gate. Incoherent operations such as cooling, preparation and detection were performed 
on the S 1/2 —> P \/2 electric-dipole transition using a 422 nm light in addition to two 
repumps at 1092 nm and 1033 nm in order to deplete the population in the D 5 / 2 and 
D 3/2 levels. State initialization to the | /S') state was done by optical pumping with 
circularly polarized light. Detection was carried out by counting fluorescence photon 
using a photomultiplier tube for a typical duration of 1 ms. State discrimination relies 
on the fact that only population in the Si/ 2 contributes to the fluorescence. From the 
statistics of the counted photons the probabilities of having zero, one and two qubits 
(for two-qubit register) in the IS 1 ) state were inferred. An example of such histogram is 
shown in Fig.lb. 

Manipulating optical qubits requires a narrow linewidth laser. As an example, the 
dephasing time of optical qubits will be limited by the inverse of the laser linewidth as 
the laser itself constitutes the local oscillator reference. Stabilized diode lasers are 
a relatively simple and cheap option for the realization of narrow linewidth lasers. 
Their use has been growing in popularity in the field of precision measurements and 
in particular, in atomic clocks where sub-Hz linewidth has been obtained [10-12], Using 
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PM fiber 



Figure 2. Schematic drawing of the narrow linewidth diode-laser at 674 nm. 

narrow linewidth lasers for coherent operations on trapped ions posses slightly different 
requirements than atomic clocks as the former is highly sensitive to the spectral purity 
of the laser at a larger bandwidth i.e. in some operations, as two qubit gate, the 
laser is tuned to the motional sidebands which are few MHz away from the carrier 
transition resonance. Therefore any spectral noise at this range will drive undesired 
carrier transitions that will reduce the fidelity of the gate. In our case, the main 
contribution to the noise at this spectral range arises from the servo loop which locks 
the master laser to a stable high-finesse cavity [13]. When the fast feedback is based on 
current modulation of the diode laser the typical unity gain, at which the induced noise 
is maximal, is at about 1-2 MHz. As we will show, a way to overcome this problem is to 
utilize the light that is transmitted trough the cavity. With a linewidth of several kHz, 
the cavity functions as a narrow optical filter which reduces the high frequency noise 
significantly. 

2. Narrow linewidth diode laser 

2.1. The master laser 

Light is generated by an external cavity diode laser (ECDL, Toptica DL100) that is 
locked to a high-finesse Fabry-Perot cavity which is made of ultra-low expansion glass 
(ULE). The ECDL, with a 10 cm long cavity, produces 10 mW of power at 674 nm 
and is enclosed in an acoustically-isolated box. The laser frequency is stabilized to a 
cavity resonance by means of the Pound-Drever-Hall (PDH) locking technique [14], A 










































Universal gate-set for trapped-ion qubits using a narrow linewidth diode laser 


5 


single mode polarization maintaining (PM) fiber guides the light from within the box to 
the optical table where sidebands at 20 MHz are added to the laser spectrum by phase 
modulation using an electro-optic modulator. Roughly 500 //W of light is sent by a PM 
fiber to the ULE cavity which is mounted on a separate breadboard. A polarizing beam 
splitter (PBS) and a A/4 retardation plate direct the reflected light from the cavity 
toward a 100 MHz bandwidth photo-detector to produce the PDH beat signal. The 
beat signal is demodulated and fed into a fast analog controller ( Toptica PDD and 
FALC modules) to optimize the transfer function of the servo loop. The output from 
the controller is fed-back into a field-effect transistor (FET) connected in parallel with 
the diode laser inducing current changes, which in turn steer the laser frequency. Due 
to the limited dynamic range of the FET, another feedback loop, implemented on a 
field programmable gate array (FPGA), keeps the average voltage of the FET fixed by 
controlling the voltage on the piezoelectric transducer that controls the ECDL grating 
angle. 

The ULE reference cavity (Advanced Thin Films Inc) has a finesse of T = 10 5 
measured using ring-down spectroscopy. The cavity is placed in a vacuum chamber 
at a pressure of 10 -7 Torr. Viton rubber pads are used to reduce acoustic vibrations. 
The free spectral range of the cavity longitudinal modes is 1.93 GHz and the cavity 
linewidth is T c = 22 kHz. The cavity is thermally isolated by the vacuum environment 
and a thermal isolation material which is wrapped around the vacuum chamber. The 
resulting peak-to-peak amplitude of the cavity frequency drift at 674 nm is about 1 
MHz over a time scale of hours and the maximal slope is about 2 kHz per minute. Since 
the laser has to be tuned to the atomic transition frequency with sub kHz accuracy, we 
scan the atomic transition every three minutes and by linear interpolation evaluate the 
cavity resonance frequency and adjust the laser frequency accordingly at the beginning 
of each experiment. 

2.2. Laser frequency and amplitude control 

The various operations on the optical qubit set a few requirements on the laser frequency 
control setup which is based on acousto-optic modulators (AOM): 

(i) Frequency detuning range up to the trap driving frequency of 21.75 MHZ 

(ii) Driving two frequencies simultaneously for the MS-gate 

(iii) Shifting all the above output frequencies to compensate for slow ULE cavity drifts 

(iv) Keeping a controlled phase relation between the various frequencies 

To satisfy all demands the use of two AOMs is unavoidable. One is in a double-pass 
configuration and the other is in a single-pass configuration. While the double-pass 
allows for a wider range of detuning it can not be used to produce a bichromatic 
beam. This is because when the beams are retro-reflected in the second pass they 
will diffract from both frequencies that drive the AOM, leading to a third beam with 
frequency shift of /j + /2 which is unwanted. Therefore, the bichromatic beams are 
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Figure 3. Schematic drawing of the various frequency sources for the single-pass and 
double-pass AOMs. The laser frequency that reaches the ions is always the sum of the 
double pass which switches between carrier and micromotion operations and the single 
pass which tunes to the carrier or vibrational sideband. 

generated with the single-pass AOM. We use a separate frequency source for almost 
each frequency in the experiment in order to control the phase relation between the 
various frequencies. The alternative of using a single frequency source and varying its 
frequency will introduce phase noise due to time jitter in the frequency updating time. 
The schematic diagram of our frequency sources is presented in Fig.3. The double pass 
AOM has three frequency sources which are all based on direct digital synthesizers 
(DDS, Novatech 409). An rf switch is used to choose between the carrier frequency and 
the micro-motion (MM) sideband frequency (used for single qubit addressing which is 
explained in the next section) separated by 21.75/2 MHz. The two frequencies are shifted 
by the base frequency electronically with a mixer where the fundamental frequencies are 
filtered by a bandpass filter. We use the base frequency to correct for the ULE drifts 
without the need to update both the carrier and MM-sideband DDSs which would 
introduce relative phase noise between them. The single pass AOM has four frequency 
sources. The red sideband (RSB) and blue sideband (BSB) generate the bichromatic 
beam for the two-qubit entangling gate. FI is used to drive carrier transition with a well 
defined phase with respect to the RSB and BSB. The last frequency source F2 is a DDS 
that updates its output frequency in real-time to drive operations with no phase relation 
requirement such as electron shelving for state detection and ground-state cooling. 

2.3. Laser high-frequency noise - servo bumps 

The bandwidth of the laser-cavity locking servo was found to be limited by the response 
of the laser diode to current modulation. Increasing the servo gain to reduce the laser 
linewidth is impeded by the growth of oscillations around the frequency of inverse 
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feedback phase. Even when the gain at this frequency is kept below unity to maintain 
stability, noise around this frequency is amplified. These amplified noise bands around 
the unity gain frequency are known also as ’servo-bumps’. In our system these servo- 
bumps are located at about 1-2 MHz and have FWHM of about 1 MHz. We observed 
servo-bumps in the spectral content of the laser experimentally in two different ways. 
The first is a direct measurement of the PDH error signal spectrum with an rf spectrum 
analyzer. This measurement result is plotted in Fig.4a. Secondly, we measured the ion 
response to an excitation of the quadrupole transition. Here, the laser frequency was 
scanned across the transition and the population in the D 5 / 2 state was measured. The 
duration of the excitation pulse was set to 100 fis, and the laser propagation direction 
was perpendicular to the axial direction of the trap to avoid coupling to the axial motion 
of the ion in the trap. The result of this measurement is presented by the blue solid 
line in Fig.4b. The fast oscillations around the central resonant peak are due to the 
100 kHz Rabi frequency which results in many oscillations during the 100/is pulse. The 
additional narrow peaks at around 0.5 MHz are due to second harmonic of a radial 
sideband of a transition between different Zeeman states. Both the Rabi spectroscopy 
and the spectrum analyzer are showing a similar feature of a 1 MHz wide response around 
±1.1 MHz from resonance. This wide response is attributed to incoherent excitation of 
the carrier transition due to the servo-bumps in laser spectrum. This spectral feature 
has undesired consequences when driving vibrational sidebands transitions that are 
typically detuned by a 1-2 MHz from the carrier. In this case the servo-bump will 
drive incoherently the carrier transition and lead to decoherence. For example, we have 
observed that the fidelity of the MS entangling gate is severely compromised by this 
effect. 

One way to eliminate the harmful effect of the servo-bump is to improve the 
performance of the servo loop by increasing its bandwidth. In this case the unity-gain 
frequency, and hence the servo-bumps, will be shifted towards a higher frequency and 
away from the motional sidebands. This can be done, for example, by introducing an 
intra-cavity EOM into the ECDL [15]. Here, We have taken another approach in which 
we optically filter those servo-sidebands from the laser spectrum [13]. Fortunately, the 
cavity itself acts as the narrow bandpass optical filter with a 22 kHz width Lorentzian 
transfer-function. Hence, the light that is transmitted through the cavity possesses 
high spectral purity with, in principle, three orders of magnitude reduction in the noise 
spectral power density around the 1 MHz servo-bumps. 

The price for relying on the transmitted light is the attenuation in the available light 
intensity. In our setup , the cavity transmission is only about 5%. While a decrease 
by a factor of two can be accounted for spatial mode-mismatch, most of the loss is 
due to intra-cavity losses. In addition, the intensity of the light that is sent to the 
cavity is limited to the mW range in order to prevent substantial thermal effects on 
the cavity mirrors (ULE substrate) that would in turn lead to thermal shifts of the 
cavity resonance frequency. Overall the transmitted light intensity is limited to few 10’s 
of 1 -iW, which is well below the required intensity for ion-qubit operations, and hence 
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Figure 4. Spectrum of the 674 nm laser as obtained by a) the PDH error signal and b) 
excitation spectroscopy on the Se->D transition with (blue) and without (red) filtering 
the laser with the cavity. The inset shows the excitation probability as function of the 
pulse duration at the peak of the servo-bump at around 1 MHz. 

needs to be amplified. To this end, we incorporate another slave diode-laser which is 
injection-locked to the filtered master-laser light. A slave diode rather than a tapered 
amplifier is chosen due to the much smaller intensity that is needed to injection lock a 
diode as compared with the intensity that is required to saturate a tapered amplifier. 

The complete master-slave setup is shown in Fig.2. The transmitted light 
polarization is restored to linear by another A/4 plate and coupled to a PM-fiber. Two 
optical isolators are placed at the input and output of the cavity to prevent optical 
interference. Between 5-10 /il'F of filtered light is injected to the slave diode laser using 
mode-matching optics. This intensity is sufficient for a stable operation of the slave 
laser with an injected current window of 0.1 m A. 

Intensity fluctuations of the 674 nm light lead to fluctuations in the Rabi frequency, 
and therefore affect any coherent operation. The origin of intensity noise is mainly due to 
fiber phase noise that is translated to intensity noise when passing through a polarizing 
beam splitter (PBS) on the way to the ion. To prevent this, the intensity of the laser 
is monitored close to the input port of the vacuum camber and stabilized by adjusting 
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the rf power that drives the AOM. 

Figure 4b compares between the excitation spectrum of the ion when excited with 
the slave laser when it is injection locked to the unfiltered master laser (blue curve) 
and to the cavity filtered light (red curve). The significant improvement in the spectral 
purity is clearly observed. The inset shows a measurement of the excitation probability 
as function of the pulse duration at the maximum of the servo-bump in the two cases. In 
the unfiltered light case already at a pulse duration of 400 jis the population is saturated 
close to 50% while in the filtered light case as little as 1 - 2 % are observed. The ratio 
between the slopes of the two curves in the linear regime is about 50 which should be 
regarded only as a lower bound due to the small signal in the filtered case which could 
be a result of power broadening together with other decoherence mechanisms such as 
magnetic field fluctuations or leakage of the 1033 nm repump light. 

2.4. Characterization of the laser linewidth 

In quantum information processing the qubit coherence time is determined by the 
stability of the relative phase between the qubit and the local oscillator. In the case 
of an optical qubit, the local oscillator that keeps track of the qubit phase, is a narrow 
linewidth laser. The physical mechanisms that cause phase fluctuations to the qubit 
and the laser are different and uncorrelated. Therefore, reaching a long coherence 
time requires that the phases of both the qubit and the oscillator are independently 
stable. While commercially available oscillators at radio and microwave frequencies can 
easily reach coherence times of many seconds (requiring fractional frequency uncertainty 
of 10 -6 — 10 -9 ), obtaining a laser with a coherence time of one second, i.e. 1 Hz 
linewidth, is a much more challenging task and requires fractional frequency uncertainty 
of 10 -14 — 10 -15 . On the qubit side, the dominant source of dephasing are fluctuations 
of the ambient magnetic field . 

One method of characterizing the spectral properties of a laser is a delayed self¬ 
heterodyne interferometer [16]. Here the laser light interferes with a delayed fraction of 
itself by using a long optical fiber to create a substantial time delay r. However, as the 
spectral resolution of this method is given by 1/r, to measure the sub-kHz frequency 
fluctuations requires hundreds of kilometers of fiber length which render this method 
impractical. The alternative is to compare against another stable reference. Examples 
for such a reference include another independent laser, an ultra stable Fabry-Perot cavity 
or a narrow atomic transition. 

Since the purpose of our laser is to drive the S\/ 2 —>■ D 5/2 transition in 88 Sr + and 
with no other tool at hand, we inferred our laser linewidth by performing a Ramsey 
spectroscopy measurements on the ion. One complication that arises in our case is 
that the relevant quadrupole transition is first order Zeeman sensitive so a conventional 
Ramsey experiment can not discriminate between the magnetic field noise and the laser 
frequency noise. Previously, custom-designed entangled state were used to differentiate 
laser phase noise from magnetic field noise [17] Here we demonstrate a custom-designed 
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Figure 5. Magnetic field noise dynamic decoupling scheme. The Zeeman shifts of two 
S 1/2 ground state sub-levels with respect to the 5D 5 / 2 ,-i/ 2 have opposite signs. Hence, 
the contributions of the magnetic field noise to decoherence in a Ramsey experiment 
can be suppressed by applying an appropriate sequence of rf 7 r pulses. 


dynamic decoupling sequence to this end. The advantage of using dynamic decoupling 
for this purpose is the simplicity of its implementation. Our magnetic field dynamic 
decoupling (MFDD) eliminates only phase fluctuations due to slow magnetic field 
noise, thus rendering the optical transition to an effective clock transition. After 
initializing to the superposition state = 1 /v / 2 (|*S'i /2 -1/2} + l-Ds/2,-1/2}) the scheme 
follows conventional dynamic decoupling methods [18] in the form of echo pulses with 
two main differences. The first is that the pulses are performed in the S 1/2 spin manifold 
and not on the optical transition as this will reverse the effect laser frequency noise as 
well as magnetic field noise. The basic idea is illustrated in Fig.5 and relies on the 
opposite sign of the Zeeman shift of the £1/2,-1/2 115/2,-1/2 transition as compared 

with the S'1/2,1/2 <->■ D5/2-1/2 transition. As long as the magnetic field does not vary 
between echo pulses, the phase that is accumulated when the ground state population 
is in the S'1/2,-1/2 due to an arbitrary Zeeman detuning, can be reversed by transferring 
the population to the Si/2,1/2 state where the phase is accumulated with the opposite 
sign. The second difference from a regular echo sequence is that the time which is spent 
in each state is not equal. This is because the magnitude of the Zeeman shift in the 
two states is different with a ratio of 1:4. The experimental sequence is illustrated in 
Fig.5b. For a Ramsey experiment of a total time T and N rf echo sequences, each 
echo sequence is of duration Ty = T/N and is composed of two rf 7 r pulses separated 
by a time T/v/5 at the middle of the sequence. Notice that in terms of the optical 
transition this scheme is identical to a conventional Ramsey sequence having only two 
7 t /2 pulse with a wait time in between where it is only the ground state part that is 
manipulated with rf pulses. We have shown, in previous work [19], that the coherence 
of ground-state superpositions can be maintained for as long as a second using several 
hundred echo pulses. In general any dynamic decouple scheme can be adapted for the 
manipulation of the ground state part. The experimental results are presented in Fig.6a. 
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Figure 6. Measurement of the laser phase correlation time by a Ramsey experiment, 
a) Fringe contrast as function of the Ramsey time with (red) and without (blue) MFDD 
sequence with N=4 as shown in Fig5b. The solid blue (red) line is a theoretical fit 
to a Bessel (Gaussian) function, b) Fringe contrast as a function of the experiment 
time in Hahn-echo sequence where a single echo n pulse is performed on the optical 
transition. Solid black line is a theoretical fit to a decaying exponent. The Inset show 
a Rabi spectroscopy scan which is consistent with the result of the Ramsey and echo 
measurements. 


We plot the Ramsey fringe contrast as a function of the total experiment time for a plain 
Ramsey sequence (blue dots) and MFDD sequence with N=4 (red dots). The contrast is 
obtained by scanning the phase of the second tt/2 optical pulse and fitting the data to a 
sine function with the amplitude and phase as fit parameters. Solid lines are theoretical 
fits to a Bessel function in the plain Ramsey case and a Gaussian in the MFDD case. 
In the regular Ramsey experiment the contrast dropped to 50% after 1.7 ms. Here, the 
50 Hz monotone magnetic field noise component dominated and the contrast followed a 
Bessel shape [20]. Applying the MFDD sequence resulted in a coherence time of 6.7 ms 
(more than a four-fold improvement) and in addition the contrast follows a Gaussian 
function shape rather than a Bessel, which implies that there is no single frequency 
component that dominated the noise spectrum. Applying MFDD sequences with N> 4 
did not extend the coherence time. Hence, we infer that at this point the laser frequency 
noise, rather than the magnetic field noise, is the limiting factor. Assuming a Gaussian 
line-shape, we infer a laser linewidth of 65(5) Hz. This result is consistent with our Rabi 
spectroscopy measurements depicted in the inset of Fig.6b. More information on the 
laser spectral content can be obtained by a standard single echo sequence on the optical 
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transition. Here, we did not apply the MFDD sequence but rather incorporated an active 
stabilization system that reduced the ambient magnetic field noise to a negligible level. 
As shown in Fig.6b. the contrast extended to 20 ms following exponential decay. Since 
the echo pulse acts as a highpass filter we infer that the coherence time measured without 
the echo is limited by low frequency noise which we assume to has a 1/f characteristic. 
If we consider the echo measurement to represent the laser fast-linewidth (white noise 
contribution) then from the exponential fit we obtained 11(1) Hz linewidth [21]. 

3. Single-qubit gate 



Figure 7. Collective Carrier Rabi oscillation on the Si/2,1/2 ** D5/ 2,3/2 quadrupole 
transition in Doppler cooled ions. From the fit (solid red line) we deduce a Rabi 
frequency of ti/2ir = 193 kHz and a single flip fidelity of 0.997(2) 


For the characterization of single qubit operations in the context of quantum 
information processing we need to distinguish between operations which are collective 
and act on all the qubits in the register, and operations that are applied to individual 
qubits. The former is similar to the case of a single-qubit register in which high fidelity 
rotation gates that go beyond current estimates for fault tolerance threshold have been 
demonstrated [22,23]. In our case, we perform single qubit rotations by driving the 
optical electric quadrupole transition. Figure 7 shows 20 consecutive collective carrier 
Rabi flops on a two-qubit register with Rabi frequency of Ut c /2n 200 kHz. The small 
reduction in the contrast can be attributed to the laser frequency and amplitude noise 
as well as the motional spectator modes as the ions are only Doppler cooled with axial 
(radials) secular frequency of 1.95(3.8) MHz. From the fitted decaying envelope we can 
deduce a single flip ( 7 r rotation) fidelity larger than 0.997(2) but a more reliable result 
would be inferred by performing randomized benchmarking [24] which is beyond the 
scope of this work. 
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In the case of more than one qubit, where single qubits need to be addressed 
individually, the experimental challenges are different. Here only a handful of 
experiments have demonstrated single qubit addressing and gates with high fidelities 
[4, 25-27]. The major challenge here is to apply a high fidelity operation on one qubit 
while preventing unwanted state changes to the neighboring qubits. The two most 
popular methods rely on tightly focused laser beams that locally induce a light shift 
on a singe ion qubit and strong magnetic gradients, both lift the degeneracy between 
the different qubits in the register and allow to individually address the desired qubit 
in the spectral domain. Our approach uses the readily available spatial variation of 
the rf electric field amplitude of the ion trap along the axial direction [28]. The basic 
principle, which is illustrated in Fig.8a., is to utilize the excess micromotion to drive 
the optical qubit transition. For the case of only two qubits, which is discussed here, we 
position the two-ion crystal such that one is exactly at the rf null (null-qubit) whereas 
the other qubit (MM-qubit), which is not at the null, undergoes micromotion at the 
trap drive frequency of 21.75 MHz. The excess micromotion that we use is along the 
axial direction of the trap and is due to the boundary conditions that are imposed by 
the endcap electrodes. In the ion frame of reference, the periodic Doppler shift due 
to excess micromotion adds sidebands to the laser spectrum. By detuning the laser 
frequency from the carrier transition by the trap rf frequency, the MM-sidebands can be 
used to resonantly drive coherent operations on MM-qubit while the Rabi frequency for 
null-qubit is largely reduced and ideally nulls altogether. The Rabi frequency for driving 
the MM-sidebands depends on the micromotion amplitude, Qmm = D c • Ji(xmm • k), 
where fl c is the carrier Rabi frequency in the absence of excess micromotion, J\ is 
the first order Bessel function, k is the laser wavevector and xmm is the micromotion 
amplitude vector. When ~x. M m is much smaller than the laser wavelength the Rabi 
frequency can be approximated by ft mm = D c /2 • xmm • k. Moreover, the reduction in 
carrier Rabi frequency due to the micromotion is second order in the amplitude and can 
be neglected. 

Fig.8b shows single qubit Rabi flops by driving the MM-sidebands. The various 
curves represent the different measured register states PO (blue), PI (black) and P2(red) 
which correspond to the probabilities of finding zero, one and two qubits in the | S) state. 
As is expected from only one qubit flopping, PI and P2 oscillate with an opposite phase 
with almost 100% amplitude while PO stays below 1% throughout the measurement. 
This means that the null-qubit remained in the |,S') while the MM-qubit oscillated 
periodically between the \S) and | D) states with a Rabi frequency of roughly 20 kHz. 
Here, the axial frequency was lowered to 1 MHz which increased the separation between 
the ions and therefore the micromotion amplitude and Rabi frequency of the MM- 
sidebands. In this measurement the single qubit flip (7r rotation) fidelity was as high as 
the collective qubits flips that are presented in Fig.7. However, the MM-sidebands gate 
is not as robust and stable as collective carrier rotations. This is due to the extreme 
sensitivity of the excess MM to the electrostatic environment where very small changes 
in the ambient electric fields change the MM-sidebands Rabi frequency and introduce 
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a) 




b) 




Figure 8. Single qubit addressing based on MM-sideband. a) Graphic representation 
of the single qubit addressing scheme. The two ions experience different micromotion 
due to their separate locations in the trap. Consequentially, the spectrum of a spatially 
homogeneous laser beam will be different in the ions frames of reference. Detuning the 
laser frequency by the trap drive frequency, results in a Rabi frequency that depends 
on the MM amplitude, b) One qubit Rabi oscillation in a two-ion crystal by driving 
the MM-sideband. The Rabi frequency is ~ 20 kHz. 


rotation errors. To mitigate the drift in the MM-sidebands Rabi frequency, composite 
pulses can be applied. As an example, the sequence X 7r / 2 TJ r X n / 2 where X, Y are the 
direction of rotation, implements a ir rotation that is first-order insensitive to small 
variation in the Rabi frequency. 

The two operations discussed above: two-qubit collective rotations and individual 
addressing of only one qubit are sufficient in the case of a two-qubit register to 
perform any single qubit gate. As an example in ref [29] we have used this scheme 
to perform complete quantum process tomography for the Molmer-Sorenscn two qubit 
gate. In principle the use of the spatially varying electric field and micromotion induced 
sidebands can be generalized to more than two ions by introducing the concept of 
dressed-state picture [28]. 
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4. M0lmer-S0rensen two-qubit interaction 

Experimentally implementing a universal two-qubit gate is typically more challenging 
than implementing single-qubit rotations. This is because two-qubit gates require the 
synthesis of bipartite interactions between qubits. Since trapped ions are separated by 
a few microns, the interaction between their qubit states is synthesized by coupling 
their internal states to their quantized collective modes of motion. Here we present our 
implementation of the Molmer-Sorensen (MS) entangling gate [30] which, together with 
a universal single qubit rotations, forms a universal gate-set for quantum computation. 
The MS scheme, besides being successful in reaching fidelities as high as 99.3% [31], 
holds some desirable properties such as resilience to finite temperature, straightforward 
generalization to multi ion-qubit gates and relative simplicity of implementation. The 
MS interaction scheme is presented in Fig.9. It consists of a bichromatic field that drives 
collective internal state atomic transitions through the red and blue sidebands of a single 
vibrational mode. While in general the MS interaction entangles the internal state of the 
ions with their motion, at certain periodic times the internal states and motion return 
to be disentangled. At those times the ions internal states undergo collective rotations 
and the motion returns to it initial state. The interaction Hamiltonian in the rotating 
wave approximation takes the following form [32] 

H in t = ^(*77e-* A Ve““* + ae iSt ){a ( y + a ( _ 2) ) + h.c.) (1) 

where O is the carrier Rabi frequency and we assume an equal intensity in each of the 
two beams, // is the Lamb-Dicke parameter, 6 = (cui — uy)/2 — v is the bichromatic 
field symmetric detuning from the motional sidebands, A = (coj, + ny)/2 — lua is the 




ss 


2(5 + v) 


Figure 9. The M0lmer-S0rensen interaction scheme. A bichromatic beam couples the 
| SS,n) to the \DD,n ) through the intermediate states \SD,n±l) and \DS,n ± 1). 
The two frequencies oj r ^ are ideally symmetrically detuned from the carrier transition 
frequency ujq by S + v, where v is the frequency of the vibrational mode through which 
the gate is performed. In addition the frequencies can be non-symmetrically detuned 
from the carrier transition by an asymmetric detuning A. 
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bichromatic field asymmetric detuning from the two-photon transition frequency, a and 
are the creation and annihilation operators of the vibrational normal mode with a 
frequency v and ct_ j+ are the ladder operators for the atomic transition. 

In the absence of asymmetric detuning, A = 0, an analytic expression for the 
time evolution of atomic populations can be obtained [33]. When starting from the 
vibrational ground state (n=0) the population dynamics takes the following form, 

P0 = 1(3 + e 2|a| - 4e |a|/2 ) + sin 2 (d)e" |a|2 (2) 

8 

Pl = \(l-e-W 2 ) (3) 

P2 = 1(3 + e 2W - 4e w/2 ) + cos 2 (#)e- | “ |2 (4) 

8 

where 

9(t) = — ^ ^(1 — sinc(hf)) and ait) = ^(e l<5t — 1) 

2 o o' 

Starting from the state | SS), a maximally entangled state |$) = (| SS) + i \DD))/y/2 
is obtained with a detuning of 8 — 2//Q and after an interaction time of T gate = 2ir/5. 
In cases where the asymmetric detuning non-zero, the dynamics can be obtained by 
numerical integration. 

4-1- Extensive 2D scans of the MS-interaction parameters 

Our experimental implementation of the MS-interaction is similar to that outlined in ref 
[33]. A bichromatic beam was produced by single-passing an AOM(AOM-S) which was 
driven by two rf frequencies (RSB and BSB DDS sources). The beam then propagated 
to the ions through a polarization maintaining single mode fiber after which it was 
focused to a waist of about 20 /xm. At the trap center each spectral component had up 
to 1 mW of power at 674 nm. At each experiment the two ions were initialized to the 
| SS) state and their stretched mode which was used for the MS-interaction, was cooled 
to the ground state (n ~ 0.05). 

As a first step, we examined the ions state population for a wide interval of the 
interaction parameters: symmetric detuning S, asymmetric detuning A and interaction 
duration t. Here the center of mass vibrational mode with frequency of u = 0.98 
MHz was used. These scans provide a visual insight to the population dynamics and its 
dependence of the various parameters of the interaction. In addition, some experimental 
uncertainties in the gate parameters can substantially effect the gate fidelity. As an 
example, uncertainties in the intensities of each component in the bichromatic beam, 
the vibrational mode frequency and light shift from other levels affect the required gate 
time and detuning. As shown in [33], the experimental calibration of all these parameters 
using the gate itself is possible. A different approach, which is more straightforward, 
however time consuming, is to scan the relevant parameters: symmetric detuning 8, 
asymmetric detuning A and time T gate and measure the two-ion response. Once the 
experimental maps are obtained the real parameters can be by extracted by comparing 
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Figure 10. Populations time dynamic for the MS-interaction as function of the gate 
detuning. For constant asymmetric detuning A = 0. Upper panels are the measured 
populations and lower panels are the calculated result of Eqs. (2)-(4). The shift in the 
y axis between the measured and calculated result is exactly the uncompensated light 
shift contribution in the measurement of the vibrational mode which was obtained by 
sideband spectroscopy. 


the measured maps to the calculated ones. Figures 10-12 present the states population 
maps, both measured and calculated, for three 2D scans: symmetric detuning vs. time, 
asymmetric detuning vs. time and asymmetric vs. symmetric detuning. 

The practical value of such maps is best demonstrated in the time and symmetric 
detuning scan which is shown in Fig. 10. The calculated populations as function of 
the gate time and symmetric detuning fully describe the measurement result up to 
a shift in the detuning of 35 kHz. This large shift is attributed to a bias in our 
vibrational frequency measurement scheme. Our sideband frequency was measured via 
sideband spectroscopy using intense beams which generated light shift on the carrier 
transition frequency which was different in the MS gate beam configuration. Therefore, 
the measured sideband frequency included a significant light shift of 35 kHz. 

4-2. Two-qubit entangling gate 

After a coarse setting of the gate parameters. Another set of quick one-dimensional scans 
were conducted in order to fine-tune the gate parameters and maximize the gate fidelity. 
First we scanned the gate symmetric detuning at the presumed gate time and set it to 
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Figure 11. Populations time dynamic for the MS-interaction as function of the 
gate asymmetric detuning A with the optimal gate symmetric detuning S/2tt = 10.5 
kHz. The calculated result shown by the lower panels was obtained by numerically 
integrating the propagator of the Hamiltonian in Eq.(l). 


a value where the populations P0~P2 ~50%. Next, the gate time was scanned with the 
obtained gate detuning. Here the gate time was set to the point where PI is minimal 
(ideally zero). The two scans were repeated several times to make sure we converge to 
the correct values. Figure 13a. presents the populations as a function of the interaction 
time. Here, we used the stretch vibrational mode . At a gate time of T gate = 140 /.(S 
the ions superposition consists of only the | SS) and | DD) states. To verify that the 
resulting state is indeed a Bell state, the off-diagonal elements of the density matrix 
must be measured as well. This was done by applying another 7 r /2 pulse with a phase <fi 
and measuring the state parity P = P0 + P2 — PI. The amplitude of parity oscillations, 
Ap, as function of <fi is equal to 2\pss,dd\- The gate fidelity for producing the Bell state 
|$) is given by i> = (P0 + P2)/2 + A p . Using the maximum likelihood method to fit 
the data in Fig.l3b, we conclude gate fidelity of F=0.985(10)11. While part (~ 5 • 10 -3 ) 
of the gate error can be attributed to imperfection in the preparation and detection 
operations, the rest is assumed to arise from error in the gate itself. Although the exact 
source of error has not been fully identified, the nature of the error has been studied 
in [29] where we have performed process tomography on five consecutive gates. The 

|| The likelihood of the measured data was calculated given a series of binomial distributions with mean 
values that follow a sinusoidal pattern with a given amplitude and phase. The amplitude and phase 
that yield the highest likelihood were used as an estimate for the gate fidelity. 
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Figure 12. Populations for the MS-interaction as function of the gate symmetric 
detuning and asymmetric detuning for a fixed interaction time t=95 //S which is the 
gate-time for the appropriate detuning, <5/27r = 10.5 kHz and A = 0. 

results were consistence with a collective depolarization channel. A possible source of 
error is off-resonance \S) | D) incoherent transition as indicated by the measurement 

in the inset of Fig.4b. 

5. Summary 

We presented a universal gate-set for qubits encoded in the electronic states of 88 Sr + ions. 
The gates were implemented using a narrow linewidth diode laser. Diode lasers have 
a relatively wide frequency noise spectrum to begin with and therefore are challenging 
when trying to reach the narrow linewidths required. In particular, we have observed 
that unlike the case of optical atomic clocks, where noise faster than the linewidth of 
the laser plays a marginal role, here, residual fast phase noise in the servo-bumps that 
overlaps the motional sideband transition can generate significant errors to quantum 
gates. We show that one can overcome this problem by using a high-finesse cavity as 
a narrow optical filter. This moderates the servo bandwidth requirement and makes 
an ECDL a compact and affordable solution as a narrow-linewidh laser for quantum 
information processing. We characterized the laser frequency uncertainty using the 
ions and found a linewidth below 100 Hz. In order to eliminate the contribution 
of the magnetic field noise from the result of a Ramsey experiment we introduced 

















Universal gate-set for trapped-ion qubits using a narrow linewidth diode laser 


20 



Gate Timers] 



Figure 13. A M0lmer-S0rensen two-qubit gate after optimizing the gate parameters, 
(a) Populations time evolution as function of the interaction time, (b) Parity oscillation 
P = P0 + P2 — PI obtained by scanning the phase of a 7r/2 pulse at the end of the gate 
at a gate-time £^=130 pS. From these result we infer gate fidelity of F=0.985(10). 


a new type of dynamic modulation sequence. This method may find use in other 
applications such as atomic clocks that are not based on a magnetic insensitive transition 
[34, 35]. We demonstrated high-fidelity collective single-qubit gates as well as individual 
qubit addressing in a two-qubit register by utilizing the inhomogeneous MM-sidebands. 
Although, this method for individual addressing is limited for only two qubits and 
not applicable for scalable quantum information processing, we still find it extremely 
useful in many cases where only two qubits are of interest. For example, we have 
successfully used this method to measure the magnetic dipole interaction between the 
spins of two trapped ions [36]. We have implemented a Mplmer-Sprensen two-qubit gate 
with fidelity of F=0.985(10), thus realizing a universal gate set for quantum computation 
with trapped ion qubits. 
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